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Abstract 

Interaction with environment may lead to the transition of quantum system from pure state to 
the mixed one. In this case, the problem of definition of entanglement may arise. In particular, 
quantitative measure of entanglement concurrence is introduced differently for pure and mixed 
states. To solve this problem, alternative definition of concurrence of mixed states is introduced in 
present work. In particular, concurrence of mixed states is determined as an average of pure states 
concurrences. Averaging is carried out over the quantum ensemble, corresponding to all possible 
realization of random variable. 

PACS numbers: 0367.Mn.89. 70. +c 
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When studying complex quantum systems, characteristics, related to the correlation of 
quantum states, often manifest themselves. This leads to necessity of introduction and con- 
sideration of entangled states. The first attempts to analyze characteristics of entangled 

n n 

states were made by Einstein m and Schrodinger [2|] in the 30-th of the 20- century. In 
spite of traditional fundamental significance of entangled state, the present splash of in- 
terest is motivated by possible application to the new physical problems. In recent years 



entangled states have become important concept of new applied disciplines such as quan- 
tum cryptography 3|, quantum theory of information {4], and physics of quantum com- 
putation [5I. It is well known, that entanglement is a specific quantum form of correla- 
;ion, which, however, possesses a number of essential differences from classical correlations 
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That is why, it describes exclusively a state of quantum me- 
chanical system and defuse classical description. Entanglement of spin systems is of special 
interest 
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15l . Il6|. This kind of problems arise in the field of Cavity Quantum Electrodynam- 
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or example. This field, in its turn, is one of the promising trends 



ics (CQED) 

for quantum computing 20, l2l|]. As is known 6|], quantum correlations are connected with 
coherent superposition states of compound quantum systems. That is why, the first simplest 
definitions of entangled states were introduced for pure quantum mechanical states |6|. In 
particular, it turned out, that quantum correlation could be described with the aid of con- 
currence. Concurrence was originally introduced as an auxiliary quantity, used to calculate 
the entanglement of formation of 2 x 2 systems. However, concurrence can also be considered 
as an independent entanglement measure 2J, |26||. The original definition of concurrence for 
bipartite two-level systems is given in terms of a special basis js]: 

|ei) = |<I>+), \e,)=z\<^-), les) =2|vl/+), |e4) = |*"), (1) 

where = (|00) ± |11))/V2 and = (|01) ± |10))/V2 are the Bell states Using 
this particular basis, the concurrence of a pure state is defined as: 



c{^) = \Y,{e^m'\- 



(2) 



The value of C(\l/) determines the existence of quantum correlation in the system . in partic- 
ular the condition C(\l/) = corresponds to the absence of correlation. The maximum value 
of concurrence is reached for Bell states and is equal to C{ej) = | ^(ei|ej)^| = 1. But pure 

i 

entangled states are only the simplest case of quantum correlated systems. Generalization 



of introduced concepts and characteristics, for the case of mixed states, is non trivial and 
often equivocal. Formally, entanglement of mixed state is determined as minimal average 
value of entanglement of ensemble of pure states, that realize the given mixed state. The 



concurrence of mixed states is given by 
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Crm.ed{p)= mf (3) 

where 

p = 5^P,|^,)(^,|, >0. (4) 



Formula ([3]) makes it possible to compute concurrence o: 



of two qubits system. In particular, as was shown in 
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■ mixed states for arbitrary states 
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Cmixed{p) = max{0, -Ai, -A2, -A3, -A4} (5) 

where Aj are the eigenvalues, in decreasing order, of the Hermitian matrix R = a/ \fpP\fp 
and p = [oy ® ay)p*{ay ® Cy). But as we have noted definition of concept of entanglement 
for mixed state is not trivial. When studying non- stationary dynamics of chaotic quantum 
mechanical systems, subjected to the action of environment, the following scenario may 
be realized: At the initial moment of time t = 0, the system being in the pure quantum 



mechanical state, during time evolution performs irreversible transition to mixed state 
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25|. In this case, is not clear which definition ([2]) or ([3]) needs to be applied. Since, at 
the different moments of time, system is in different states. In this work we shall try 
to introduce alternative, more adequate for non-stationary chaotic systems, definition of 
concurrence. With that end in view we shall consider generalized Jaynes-Cummings model. 



introduced in the work (26| . This model describes a system of two interacting qubits, under 
the action of cavity field. Mixed state is formed due to the averaging over quantum ensemble 
corresponding to all possible realizations of stochastic parameter. As it will be shown below, 
random trajectory of a system plays a role of this parameter. It should also be noted, that the 
values of mixed state parameters Pi in (j4]) are determined exactly in our case. That is why 
we can not consider them as free parameters, with respect to which the expression ([3]) should 
be minimized (As requires common definition of concurrence of mixed states 6(] ) . This is 
one more confirmation of necessity of introduction of alternative definition of concurrence for 
non- stationary chaotic systems. As we have noted, the system we are interested in, presents 
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itself the generalized Jaynes-Cummings model for the case of two interacting two-level atoms 



26l | . The Hamiltonian of the problem under study has the following form: 



p2 ^ ^ 

H = \- Hs + HsB + Hb, 

2m 

Hs = MS! + SI) + ^{StS^ + S^S+), (6) 
HsB = -9ix){iSt + St)b + iS{ + 

where h = 1, is the Zeeman frequency of the spins located in the field inside the resonator, 
f2 is a constant of interaction between spins in frequency units, 3^2^ S^2 spin operators 
of two- level atoms. The last term in iQ is the Hamiltonian of the field: Hb = ujfb~^b, 6"*" and 
b are respectively the photon creation operator and the annihilation operator. As was shown 



m 



26l ] , due to nonlinear dependence of the interaction constant on the system's coordinate 



261 . As a result, the problem 



g{x) = gocos{kx), dynamics has very complex and chaotic character. Time dependence 
x{t) itself may be considered as random Wiener process. This makes a definite infiuence 
upon quantum mechanical time evolution of the system 
is reduced to the stochastic Schrodinger equation, having the following form in interaction 
representation: 

.=^^^l*W>, (7) 

where is a potential, stochasticity of which is provided by chaotic time dependence of the 
coordinate of the system x{t) : 

V = QiS+S^ + + ujfb+b - go cos{kfx) {{S+ + S^)b + {S{ + (8) 

Assume, that at the time zero {t = 0), the wave function represents itself the direct product 
of system and fields wave functions: 

= 0)) = system) ® \^ f^eld) (9) 



Here 



{^system) = Coo|00) + Coi|01) + Cio|10) + Cn|ll), 

l^/^eid) = 5^Vr„|n). (10) 
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Due to interaction ([8]), only following transitions between states are possible: 



|0,0,n + 1) ^ |0, l,n), |0,0,n+ 1) ^ |l,0,n) 
|0,l,n) ^ |l,l,n-l), |l,0,n) ^ \lA,n- 1) 
|l,0,n) ^ \0,l,n). 



(11) 



On the basis of equations (fTTl) . we shall search for the solution of equation ([7]) in the following 
form: 

I^W) = XlC'oAn+i|0,0,n+l) + ^C'o,o,„|0,l,n)+^C'i,o,n|l,0,n)+^C'i,i,„_i|l,l,n-l). 

n n n n 

(12) 

Substituting Eq. f|T2l) into Eq.Q and taking into account Eq.®, one can deduce system of 
differential equations, for time dependent coefficients of resolution of wave function. Cor- 
responding calculations takes up a lot of space and for brevity we don't bring them here 
|26|. But here we note that, time dependence of the coefficients of wave function in (fT2l) are 
determined by the functional: 

t 

i f ui(t')dt' 

Q[uj{t)] = e (13) 

where 

(14) 



u}{t) = ^/2{2nTT)goCos{kfx{t)). 



On the other hand, as was shown in j25t], because of the chaotic dynamics, x{t) may be 
treated as a classical random process. In this case, for the determination of the systems 
state it is necessary to average the functional in f|T3|) by all realizations of the stochastic 

Tom pure state to the mixed one. 



variable x{t). This leads to the transition of the system 



Since non-diagonal matrix elements of density matrix 25 



Pijkl P system ® P field 

after statistical averaging, decays in time 
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(p) ~ (e2^™Q"i[^(t)]) = e^^^^exp ( - 1 [^Erf{t^o)) 

2\ ao 



(15) 
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where Erf{. . .) is an error function 27|. Thus the time of transition from pure to mixed 
state is completely determined by the width of autocorrelation function of random variable 
xit): 



{x{t + T)x{t)) 



(17) 



In the time interval < t < ^faofTx^ the state of the system may be considered as pure. The 
transition to mixed state takes place at t > \/ olqI'k 26 |. In order to define level populations 
of the spin subsystem, we summarize all possible contributions from the different field states 

m 

oo oo 

W{t, |11|)) = W{t, |10|)) = ^ |Ci,o,„-i(t)P, 

n=0 n=0 

oo oo 

W{t, |01|)) = |Co,i,„_i(t)P, W{t, |00|)) = Y l^o,o,n+iWI', (18) 

n=0 n=0 

Using Eq. (|T8l) . the wave function of spin subsystem may be presented in the following form: 



\'^system) = (t, + (t, \01))\01) + ^/W (t, \10))\10) + ^/^^(t, 1 00) ) 1 00) (19) 

Not, that quantities W{t, \i,j)) contain terms with random time dependence. But in time 
interval < t < -^/ao/vr, the system is still in the pure state 26[. That is why, for determi- 
nation of entanglement of spin system definition for pure states ([2]) should be used. Taking 
into account (flQl) . we obtain: 

C{\^{t))s) = 2WWit, \ll))Wit, |00)) - ./Wit, m)Wit, |10))| (20) 

After lapse of time t > ^./oq/tv, random character of quantity x{t) becomes of great 
importance. Therefore ( l20l) should be averaged over all possible realizations of random 
variable. 

In order to do this, one should replace the quantities W(t, by their statistically 

averaged values {W{t,\i, j))) . Formation of mixed state in the system, is manifested in 
zeroing of time dependent interference terms in {W{t, Since this means also zeroing 

of non-diagonal matrix elements of the density matrix (16). As a result, statistically averaged 
concurrence of pure state C(|\I^(t))) will correspond to the concurrence of mixed state: 

Cmixedil'^it > vW^)) = {Cpurem)) (2l) 

The expression f l21l) may be considered as alternative definition of concurrence of mixed 



states. As a result, after decay of interference terms t » ^Jlx/J^\, we obtain: 

Crmxed = ^^WnWoo - VWo,W,o\, (22) 

where 

W,, = {W{t > vW^lu)) 



. Quantities Wij are not arbitrary constants, by which the expression (l22l) needs to be 
minimized. This is principal difference of our definition (21) from the standard one (3). 
Standard definition has nothing to do with the mechanisms of formation of mixed state and 
imphes all possible realizations of it. While, in our case mixed state is exactly defined by 
the features of the dynamic. Explicit expressions for Wij were obtained in 26|. It turned 
out that quantities Wij are expressed in terms of generalized hypergeometric functions and 
depend on average number of photons inside of cavity Wij{n). 
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2 I 1 ^ exp{—fi)y^Erfi{n] 
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3075 [2v^(-16n + exp(-n)(3 - 14n - 16^^)) - 3^/^r{l - 2nfErfi{Vfi)] 
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(25) 



Explicit form of dependence Wij{n) gives us possibility to start from essentially quantum 
domain fi ~ 1 and see what happens in the semiclassical limit Wij{n). In other words, let 
us proof, how realistic the situation is. 



In quasi-classical limit n >> 1, the form of coefficients is oversimplified 26 1. 
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W{t > y/ao/n\ll)) = W{t > vW^lOO)) = W{t > vW^|10)) = W{t > vW^|01))- 

(26) 

As a consequence Cmixedin » 1) = 0. The absence of quantum correlations in quasi- 
classical limit proves correctness of definition (21). 
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